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Abstract:  The  probability  distributions  of  the  amplitude  of  the 
envelope  energy  spectrum  of  both  Gaussirn  noise  and 
amplitude  modulated  Gaussian  noise  are  derived.  From 
these  probability  distributions  the  probability  of 
detection  of  the  modulation  is  calculated  by  comparing 
the  distribution  at  the  modulating  frequency  to  the 
distribution  at  an  adjacent  frequency.  Confidence 
percentages  as  a  function  of  modulation  factor  are 
also  calculated  and  compared  with  experimental  results. 


UNCLASSIFIED 


lit 


TABLE  OF  CONTENTS 


Page 

Acknowledgments  .  .  ,  .  .  .................  li 

List  of  Figures  .......................  iv 

List  of  Tables  .......................  v 

I.  INTRODUCTION  ....................  i 

II  SPECTRAL  DENSITY  FOR  A  SQUARE  LAW  DETECTOR  IN 

RESPONSE  TO  NARROW-BAND  GAUSSIAN  NOISE  ......  3 

III.  THEORETICAL  DISTRIBUTIONS  ...............  9 

3.1  Unmodulated  Gaussian  Noise  ...........  9 

3.2  Amplitude  Modulated  Gaussian  Noise  .......  16 

IV.  EXPERIMENTAL  ANALYSIS  .................  22 

V.  SUMMARY  ........................  39 

51  Results  and  Conclusions  .............  39 

52  Areas  for  Further  Study  ............  40 


BIBLIOGRAPHY 


41 


iv 


LIST  OF  FIGURES 


Figure 

Page 

2 . 1 

Square  Law  Detector  . 

...  4 

2.2 

Spectral  Density  of  Gaussian  Noise  ... 

6 

2.3 

Spectral  Density  for  a  Square  Law  Detector 
Given  a  Narrow-Band  oaussian  Noise  Input 

7 

4.1 

Bl^ck  Diagram  of  Experimental  Setup  ..... 

...  23 

4.2 

X-Y  Recording  of  the  Analyzer  Output  .... 

...  25 

4,3 

Analyzer  Output  Distribution  for  m  *  0  .... 

...  27 

4.4 

Analyzer  Output  Distribution  for  m  ■  0.025 

.  .  .  28 

4.5 

Analyzer  Output  Distribution  for  m  •  005  .  , 

.  .  29 

4,6 

Analyzer  Output  Distribution  for  m  *  0.075 

30 

4.7 

Analyzer  Output  Distribution  for  m  *  0,1  . 

...  31 

4.8 

Confidence  Percentages  for  t  ■  0.28  volts  . 

...  34 

4,9 

Confidence  Percentages  for  t  *  0,22  volts  .  . 

...  35 

4.10 

Confidence  Percentages  for  t  *  0,20  volts  .  , 

,  .  36 

4  11 

Confidence  Percentages  for  t  «■  0.16  volts 

...  37 

V 


Table 

4.1 


LIST  OF  TABLES 


Page 


Probability  of  the  Signal  Being  Above  the 

Threshold  .......  33 


CHAPTER  I 


INTRODUCTION 

In  recent  years,  the  study  of  the  characteristics  of  complex 
noise  signals  has  become  increasingly  mere  important  in  the  fields 
of  communications  and  acoustics.  Several  mathematical  surveys  of 
this  subject  have  been  written  by  Rice  [1]  ,  Davenport  and  Root  [2], 
and  others  Although  these  papers  cover  many  ui  me  properties  ui 
noise,  none  is  involved  with  amplitude  modulation  of  Gaussian  noise 
which  is  an  important  consideration  m  the  area  cf  signs.  detection. 

The  surveys  mentioned  aoove  contain  the  basic  material  needed 
for  this  study.  Rice's  definitions  of  narrow  band  Gaussian  noise 
and  the  envelope  of  narrow* band  jaussian  noise  are  used  throughout 
this  thesis  Of  special  interest  is  the  treatment  by  Doth  Rice  and 
Davenport  and  Root  of  the  output  of  square  law  detectors  If  the 
input  to  such  e  device  is  narrow* band  Gaussian  noise,  then  the 
output  la  proportional  to  the  square  of  the  envelope  The  authors 
treatment  11,2)  Includes  deriving  the  expression  for  the  spectral 
density  output  of  the  detector  for  a  Gaussian  noise  input. 

In  this  study  the  frequency  spectrum  distribution  of  the 
output  cf  a  square- lev  detector  will  be  calculated  This 
distribution  can  then  ba  related  to  the  power  spectral  density 
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that  has  been  derived  in  the  literature  Then  the  probability 
density  of  the  frequency  spectrum  will  be  derived  for  an  input 
of  amplitude  modulated  Gaussian  noise  to  the  square  law  detector. 

For  this  case,  a  spike  will  appear  at  the  modulating  frequency  of 
the  spectrum  and  the  amplitude  distribution  of  this  spike  will  not 
be  the  same  as  the  distribution  of  the  spectrum  at  adjacent 
frequencies.  The  amplitude  distribution  of  the  spike  will  be 
derived  separately  and  will  be  ot  primary  importance  in  finding  the 
probability  of  detection  of  the  modulation 

The  frequency  distribution,  once  derived  wil*.  be  modified  to 
a  special  case  so  that  the  theoretical  results  can  be  compared  to 
the  results  found  by  experiment  The  modification  is  a  function  of 
the  spectrum  analyzer  used  for  the  experiment  For  the  experiments 
done  in  this  study,  a  Quan  Tech  Wave  Ansiysc:  Model  30*+M  will  be 
used  Th*  theoretical  results  wi.l  be  modified  so  that  they  can  be 
compared  with  the  output  of  this  particular  analyzer 

From  the  final  frequency  distribution,  it  wil»  be  possible  to 
establish  a  threshold  level  at  the  modulating  frequency  such  that 
the  probability  of  the  output  being  above  that  threshold  for 
unmodulated  Gaussian  noise  Is  a  set  Using  this  threshold, 

the  probability  cf  the  output  being  above  that  threshold  for 
modulated  noise  csn  be  found  and  plotted  as  a  function  of  the 
percentage  of  modulation  Tt  it  these  plots  for  different 
threshold  levels  that  vil,  be  the  mein  result  ot  this  research 


CHAPTER  II 


SPECTRAL  DENSITY  FOR  A  SQUARE  LAW  DETECTOR  IN  RE3FONSE 
TO  NARROW- BAND  GAUSSIAN  NOISE 

Several  articles  have  been  written  by  various  people, 

Including  Rice  and  Davenport  and  Root,  in  which  an  expression  has 
been  derived  for  the  power  spectrum  of  the  square  of  the  envelope 
of  Gaussian  noise.  This  is  done  oy  finding  the  power  spectrum  of 
the  output  of  a  system  shown  in  Figure  2  consisting  of  a  square 
law  device  followed  by  a  low  pass  filter  If  the  input  is  narrow¬ 
band  Gaussian  noise;  nit),  then  according  to  Davenport  and  Root,  n.t) 
can  be  expressed  as: 

n't)  «  Eit»  cos  *  tf't  J  ,  il 

where  Eit)  is  the  envelope  of  the  Gaussian  ncise,  tr  *•  wjln  is 
the  center  frequency  of  the  input  spectral  density  and  0  t  is  the 
input  phase  The  output  of  the  square  law  device  wcuid  be, 

2  .  2 

Slt>  *  *  jf cos  [2fe>  t  *  49  t'i  2> 

ll  c 

If  the  bandwidth  of  the  noise  is  narrow  compared  with  the  center 
frequency,  the  spectral  density  of  the  two  terms  in  Equation  2 
will  not  overlap  Then  if  Sit,  is  passed  through  an  idee,  low 
pees  filter,  the  output  of  that  filter,  sit),  will  be 

.2 

.  1..L1. 
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s(t 


(3) 


This  can  be  normalized  by  making  a  *  2  so  that  the  output  cf  the 
filter  will  be  equal  to  the  square  of  the  envelope.. 

Assume  the  spectral  density  of  the  narrow  hand  Saussian  noise 
to  be  as  shown  in  Figure  2  2,  Then  Davenport  and  Root  ( 2 j  nave 
found  the  spectral  density  of  the  output  of  the  square  law  detector 
This  Is  i''ven  by  the  following; 

S  *  4a2N  V  &  .i  ♦  4a2N  2W' 1  -  f  i/W»  4 

z  o  c-  I 

where  N  is  the  height  of  the  spectral  density  of  the  input.  W  is 

c 

the  bandwidth  of  the  input,  and  &>f)  is  the  Dirac  delta  function 
is  shown  in  Figure  2.3  as  a  function  of  frequency, 

It  will  be  convenient  in  this  study  to  further  normalize  the 
amplitude  of  E  t*  so  that  the  height  cf  the  continuous  part  of  S 

z 

at  f  *  0  is  one.  This  will  nGt  change  the  results  and  wii:  help 
si7.pl if y  the  equations  derived.  The  constant  factor  would  drop  out 
when  comparing  the  amplitude  of  adjacent  frequencies  of  the 
frequency  spectrum  To  normalizes  let 

4a  2N  2W  *  1  «  5  " 

o 

With  this  normalizing  factor  Equation  4  can  be  rewritten  as  the 
following 


S  -  W  51  f  -  +  ,1  -  |  f  j  /V  i  1 6; 

zn 

When  the  input  noise  is  modulated  by  a  sine  waves  the 


expression  for  rhe  modulated  noise  would  be  of  the  form 
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•,  i  -i  ro  cos  w  t'j  n(t)  where  m  is  the  modulation  factor  and  ou  ia  the 
m  m 

modulation  frequency.  It  is  then  correct  to  say  that  the  envelope 

of  the  noise  13  modulated  in  the  same  way.  Therefore,  the  square  of 

2  2 .  . 

the  envelope  of  the  modulated  noise  is  (1  +  m  cos  w  t  E  it' 

ia 

The  type  of  modulated  noise  that  is  of  primary  interest  in  this 
study  is  that  of  wide  band  modulated  noise  that  is  then  passed 
through  a  narrow-band  filter.  The  modulation  produces  sidebands  in 
the  power  spectrum  which  will  have  impulses  at  plus  and  minus  the 
modulating  frequency.  The  problem  is  then  to  find  the  probability 
of  the  amplitude  of  this  impulse  being  above  a  set  threshold  as 
a  function  of  modulation  factor. 


CHAPTER  III 


THEORETICAL  DISTRIBUTIONS 


If  the  Gaussian  noise  under  consideration  is  cf  finite 
duration  of  time,  then  the  amplitude  of  the  energy  spectrum  of  the 
envelope  squared  will  be  a  random  value  with  a  certain  distribution 
If  this  amplitude  distribution  is  found,  then  its  variance  can  be 
related  to  the  average  power  of  the  envelope  squared  as  found  in 
Chapter  II.  Then  using  a  modified  distribution  which  assumes 
amplitude  modulated  noise  as  an  input,  detection  probabilities  as 
a  function  of  modulation  factor  can  be  calculated 

To  begin,  narrow-band  Gaussian  noise,  n , t > ,  has  been 
represented  by  Rice  in  the  following  way 


ni  t)  ■  a(t  »  cos  t 
c 


bi  t)  sin  u>.t 


where  a.t;  and  b(t;  are  both  i-aussian  functions  of  time  and 

f  »  u,  /2ir  is  the  midband  frequency  of  the  noise  The  square  of 
c  c 

the  envelope  is  then  defined  as  the  following: 


E2U> 


b2*  t) 
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Noise  of  bandwidth  W  can  be  specified  by  taking  samples  at 
every  ^  seconds.  If  e  record  of  noise  T  accords  in  duration  is 
considered,  then  2WT  degrees  of  freedom  are  obtained  by  casing 


10 


the  samples  at  intervals  of  ^  seconds  In  representing  the  noise 

as  m  Equation  1 ,  WI  samples  of  a(t>  and  WT  samples  of  b< t)  will 

also  specify  the  noise  and  give  2WT  degrees  of  freedom.  The 

amplitude  of  the  samples  of  a.'t)  and  bs  t )  will  have  a  Gaussian 

distribution  with  zero  mean  and  with  variance  ;  .  If  the  power 

of  the  noise  m  its  passband  is  N  ,  then  a  1  *  2N  W . 

r  etc. 

The  next  step  is  to  convert  the  samples  o*  alt'*  into  the 
frequency  domain,  Let  a.t)  be  represerted  by  its  samples  a*,  in 
Equation  9, 


ait  i 


sin  n  Wr  -  WT-r  1 

3wt  i.wt  -  wr*i 


l''9' 

where  a,,  a.;  ,  a_  are  the  magnitudes  of  the  WT  samp.es  of 

i  Z  WT 

a(t>.  If  Mi)  represents  tne  Fourier  transform  of  a  t-,  then  toe 
Fourier  transform  of  both  sides  of  Equation  9  can  be  taken  to  obtain 
the  following 


A  f>  T  ^  (a  *  a,  ecs  irf/W  la.  sen  2* f/W 
W  i  i  *. 


*  a,  ecs  Ir-t/W  •  ja.  sin  **nf/W  • 

10 

Tne  samples  of  Atf,  ate  located  at  interval*  i/T  ever  the  band 
-W/2  <  f  <  W/2  There  are  WT  aaopiea  >c  that  they  located  at 
f  »  n/I  where  n  »  *  ^  ^  1  ■  .  ^  The  value  of 

A* |p  ‘  from  Equation  10  i  * 


li 


A($) 


WT 

i  r 

W  L 

k*l 


a.  cos 

k 


(k-n2nn 

WT 


wr 


a,  sin 

k 


Ul  l2lx& 


WT 


i  l 


k*I 


The  distribution  of  a  sum  of  n  Gaussian  variables  with  zero  mean 
2 

and  variance  o  is  Itself  Gaussian  with  zero  mean  and  variance 
2 

equal  to  no^  .  Therefore L  the  samples  of  Ai~  have  an  amplitude 

2 

distribution  which  is  Gaussian..  If  represents  the  variance  of 

the  real  part  of  A(^)  and  represents  the  variance  of  the 

imaginary  part  of  A(^),  then 


WT 


cos 


2  jk^jjn 

WT 


k-1 


and 


i. 

2 


cos 


WT 


1/ 


l 


k-1 


ias 

WT 


k 

2 


CO*  ik'l 


iaa] 

UT  I 

"♦  I 


K  T 
0 


<  1? 


Nov  let  an  represent  the  samples  of  A  ^ 
samples  that  would  be  obtained  In  converting  bit 


•nd  ft  represect  th: 
n 

>  into  the  trequency 


domain.  Then 


a  ’  A  cos  in  t  *  1  A  sin  a)  t 
n  n  n  J  n  n 


4 14  * 


B  *  B  ccs  <»j  t  ■*  j  B  sin  u  t 
n  n  n  n  n 


%5 


where  A  A  *  B  „  5  are  Gaussian  variables  with  zero  mean  and 
a  n  n  n 

variance  N^T  as  found  in  Equations  i2  and  11  To  rind  the 
distribution  of  the  envelope  squared  m  the  frequency  domain,  the 
distribution  of  Aff)  convolved  with  itself  plus  B  f.  convolved  w«.th 
itself  must  be  found.  Let  F  f  be  the  result  of  the  sum  of  the 
convolutions.  Then  at  a  particular  frequency  n/I, 


WT/2 

*  y 


ia  a  *  p  p  I 

m  m-n  m  m-n 


nyo  , 


m 


WT 


where  *  signifies  the  complex  conjugate  ot  the  variable. 
Then  using  Equations  1^  and  15- 
WT/2 


*  )  ."‘A  cos  t  *  i  A  sm  t  A  cos  *  t- 

v  j  m  m  m  mmn  mn 


m 


ML 


jA  sin  w  t 
J  m-n  m  n 


*•  <6  cos  w  t  *  iB  sin  *  t»tB  cos  *  t  *  jft  sin  u,  t 
m  m  m  mmn  m  n  m  n  m-n 


m  n 


C~~* '  1  8  8 

)  [s'A  A  +  A  A  +  B  i  t  B  B 
/  ,  m  m-n  m  m-n  m  m  n  m  m-n 


COS  uj  ‘  tu  1 1 

m  m  n 


8  8  3  8 

+  (A  A  -A  A  +  B  B  -  B  B  1  cos  (u  -t  u>  t 
oi  m-n  m  m-n  oi  m-n  in  m-n  in  m*n 


A  A  ♦  i  B 


c 

1 

8 

• 

m  m-n 

m 

* 

t 

t'A  A 

-  A  A  ♦ 

3  B 

•s 

a 

a 

m  m-n 

m 

B  B  )  sin  iu*  *  u  t 
m  m_n’  m  m  n 


m-n  m  m  n 


m  r  n 


nfO  .  ' i 7) 

Each  product  of  variables  has  a  distribution  of  the  product  of  two 

Gauaslans  with  zero  mean  and  a  variance 

4 

0  2  -  — !-  •  H  'l2  .  .6 

"  4W2 

This  is  derived  from  the  general  theorem  that  states  if  twc 
Independent  variables  have  means,  and  y.  and  variance, 

and  o respectively,  then  the  distribution  of  the  preduct  cf  the 

.  .  .  .  1  2  2  2  2  2 
two  variables  has  mean  y  a,  and  variance  j,  o.  *  u  o,  ♦  y 

i  2  .  i  1  i  2  1 

In  this  case  u  *  y.  *  0  ana  *  ^  *  N  T 

1  2  1  2  O 

Equation  17  shows  that  f  consists  of  a  sum  of  i6"'iW  •*  ;  n/T 

bQ2  T* 

samples  Sil  with  zero  mean  and  variance  -2-——*  .  By  assuring 
ourselves  that  ;6I(W  -  in/Tj)  la  a  large  number,  the  central  limit 
theorem  can  be  used  to  apprexiaate  the  distribution  cf  Kjf  .  The 
only  reatriction  in  suking  16T:W  *  | n/T ! )  large  enough  is  to  mast 


14 


sure  that  the  frequency  of  interest  is  near  zero  as  compared  with 


W  or  *  W.  The  distribution  of  F »’“>  can  be  approximated  by  a 

2 

Gaussian  distribution  with  zero  mean  and  a  variance  o  r  where 

n 


2  2 
N  T 


o  *  16T/W  *  I  n/T  |  ) 
n  1  4 


2  3 

4  N  T  cW 
o 


-4N2T3Wa  •  )  £  1  /W)  „  «'  19  > 

o  1 

Since  a  finite  duration  of  time,  T,  is  being  dealt  wl'h 

2 

o  is  then  equal  to  the  average  energy  of  the  envelope  squared  in 
n 

the  bandwidth  1/T.  To  normalize  this  variance,  the  energy  must  be 

related  to  the  power  spectrum  shown  In  Figure  3  From  Equation  6 

the  average  power  per  cycle  at  a  frequency  f^  is  1  -  f  /W  -  The 

power  over  1/T  cycles  wll*  then  be  ;i  *  f^/W*  j  Therefote.  the 

energy  in  a  bandwidth  1/T  around  the  frequency  f^  la  (\  •  f^/W  ^ 

2 

Now  Equation  19  can  be  normalized  Oy  setting  or  at  f^  equal  to 
Vl  -  f^/W;. 


4N  2  T3  W  *  1  •  f  /«;  *  a  -  f  /W; 
0  XX 


Zero  frequency  is  •  epeciai  caae  end  must  be  treated 
separately.  I# writing  Equation  lb  for  n  •  0  give*  us 


i  5 


WT/2 


F(0)  *  )  {a  a  *  +  P  8  *] 

/  ,  mm  mm 


C2i> 


.  EL 


Again  using  Equations  14  and  15: 


WT/2 

V  •  ' 

F(0)  “  )  [(A  cos  u  t  +  jA  sin  u  tKA  cos  w  t  *  JA 

/  t  is  n  in  mm  in  in 


at 

■  "  -  2 


sin  u  t ) 


+  (B  cos  w  t  +  IB  sin  u  t)(B  cos  u  t  *  IB  sin  o»  t)  j 
n  o  m  tn  m  m  tnro 


WT/2 

1  I 


2  '2  2  ’2  2  "2  2  '2 
„  ,  [A  +  A  +  B  +  B  +(A  -A  +B  *B  > 

2  /  .  1  a  m  m  m  m  m  m  m 


XL 


cos  2b)  t] 

n 


22) 

Since  each  variable  is  Gaussian,  the  distribution  of  the  square  of 

the  variable  is  a  chi -squared  distribution  with  mean  equal  to  N^I 
2  2 

and  variance  2N  T  .  Again  suaning  over  these  variables  and  applying 
o 

the  central  limit  theorem,  me  can  approximate  the  distribution  of 
FtO)  as  a  Gaussian  distribution  with  a  mean 

t 

°fr  2 

u  -  4WT  Hr*  -  2  II  T  W  23, 

2  o 


8VT  i  2 


,  m 

4  # 


4T3  N  2W 
o 


v2a> 


and  variance 


16 


This  distribution  can  now  be  normalized  with  the  heio  of 
Fquation  20,  to  obtain 


and 


4T3W 

at3w 


1 


V25> 


i  26 ,. 


Now  the  amplitude  distribution  of  the  envelope  squared  in  the 
frequency  domain  Is  specified  at  any  frequency.  The  distribution 
Is  Gaussian  at  all  frequencies  with  a  mean  of  zero  and  variance  of 
(1  -  I  f|/W)  at  any  non-sero  frequency.  At  zero,  the  mean  is ~\J\H 
and  the  variance  Is  1. 


3.2  AapHw4i  N9i« 

Now  the  distribution  of  the  square  of  the  envelope  must  be 
found  given  that  the  input  noise  is  amplitude  modulated  by  a  sine 

I 

wave  of  frequency  f  ■  The  nolae  n  (t)  can  be  represented  in  the 

B 

following  way, 

« 

n  (t)  •  nit.  tl  +  m  cos  v27; 

Substituting  for  n(t>  from  Equation  7, 

♦ 

n  (t)  •  a(t)  (l  ♦  m  cos  w  t)  cos  w.c  -  b(t)  tl  -t-  »  cos  w  tj  sin  w.t. 

a  1  mi 

128) 
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t  i 

Then  from  Equation  8,  If  E  (t)  is  the  envelope  of  n  it), 


[E*(t)r 


[a(t)  (1  4  TT  cos 


4-  f b(  t )  (1  4  m  cos 


u)  t)) 
m 


(29/ 


Converting  a(t)  (1  +  m  cos  w  t)  Into  the  frequency  domain  will  give 

m  " 

A(f)  as  In  Equation  10  plus  two  sidebands  of  magnitude  A(f)  m/2. 

That  is,  if  the  samples  of  A(f)  have  zero  mean  and  variance  NqT  as 

In  Equations  12  and  13,  the  samples  of  the  sidebands  will  have 

2 

zero  mean  and  variance  N  m  / 4. 

o 

Now  let  the  samples  of  A(f)  be  denoted  by  a  and  the  samples 

of  the  upper  sideband  by  a  and  the  samples  of  the  lower  sideband 

2 

by  a  To  find  [a(t)  (1  4  m  cos  uj  t)]  tn  the  frequency  domain, 
nZ  m 

we  must  convolve  the  frequency  spectrum  of  a(t)  (1  i-  m  cos  ^  t) 

tn 

with  itself.  If  G(“r)  is  the  result  of  this  convolution,  we  will 
have: 


WT/2 

Y"1  *  *  * 

)  [a  a  4  a  ,a  .  4  a  »a  , 

/  ,  m  m-n  ml  ml-n  m2  mZ-n 

WT 

"  “  2  ’ 


tq  ,a  4  a  a  .  4  a  ,a  4  a  a  » 

ml  m-n  m  ml-n  m2  m-n  m  m2-n 


amlam2-n 


am2aml-n^ 


(30) 

Proceeding  similarly  to  the  case  with  zero  modulation,  G(;j“)  consists 
of  72T(W  -  | n/T | )  samples  with  zero  mean.  Applying  the  central 
limit  theorem  to  the  sum,  G(:p)  will  be  Gaussian  with  zero  mean  and 


variance 
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N  2T2  2N  2lV*  4N  2  2  2  2N  2T2m4 
8T(W-  I  n  /  T  | )  -V+"V-  +-V“ 


2 

•  2No2T3W(1  -  I f|/W)(l  f  j-)2 


Since  the  same  procedure  can  be  used  for  the  samples  of  b(t),  the 
distribution  of  a  sample  of  the  square  of  the  envelope  in  the 

I 

frequency  domain,  F  (f),  will  be  Gaussian  with  zero  mean  and 

variance  a  2  equal  to; 
m 

2 

a  2  =  4N  2T3W(1  -  |f|/W)(l  +  ~)2  .  (.32 

m  o  1  2 


Now  we  must  take  into  consideration  that  the  convolution  of 

the  square  of  the  envelope  will  give  spikes  at  zero  frequency, 

at  the  modulating  frequency,  and  at  twice  the  modulating  frequency. 

The  spike  at  zero  arises  when  A(f)  is  convolved  with  itself  and 

when  either  sideband  is  convolved  with  itself.  The  spike  at  f 

r  m 

arises  when  A(f)  is  convolved  with  either  sideband  and  the  spike  at 

2f  arises  when  one  sideband  is  convolved  with  the  other, 
m 

First,  consider  the  spike  at  zero  frequency. 


it  i?  it 

[a  a  t  a  ,  fa  , 
mm  ml  ml  m2  m2 


it  'it  it  it  it  it 

+  a  ,a  +  a  a  ,  f  a  ,a  +  a  a  +  a  ,a  -  +  a  ,,a  ,  ] 

ml  m  m  ml  m2  m  m  m2  ml  m2  m2  ml 


The  first  three  terms  give  distributions  which  are  chi- squared,  while 


the  remaining  terms  have  a  distribution  of  the  product  of  two 
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Gaussian  variables.  Summing  over  these  variables  and  applying  the 
central  limit  theorem,  G(0)  can  be  approximated  by  a  Gaussian 
distribution  with  mean 


|i 


»  2WT 
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'  N  T 
2 


+  2N 
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N  T2W 
o 


(1  +  f") 


<34) 


and  variance 
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64 


-  2nVv  {1  +  f-]2  <35; 

If  we  now  take  into  account  the  terms  for  b(t),  both  the  mean  and 

I 

variance  of  F  (0)  will  be  twice  what  was  calculated  in  Equation  34 

» 

and  35.  Therefore,  the  distribution  of  F  (0)  is  Gaussian  with  mean 


u  -  2N  T2W  (1  +  m2/2)  1 36) 

om  o 

and  variance 

2  2  3  2  2 

a  -  4N  Vw  U  +  m  /2 y  137, 

om  o 


For  the  spike  at  f  •  f  ,  Equation  30  can  be  written  as 

in 


c<f  ) 

in 


WT/2 

)  la.Vl  I  *  I  *  1 

—  »  m  m 

.  .  HI 

2 


*  *  *  *  * 

•t-or.a  f.+aa.-.  +  a.a  a  o  ,  .  .  + 

ml  m*f  T  m  ml-f  T  m2  ■- f  T  m  m2 - f  T  ml  m2- f  T 


a«2a.l-f  T3 


38) 


The  fifth  and  sixth  terms  will  have  a  chi-squared 


distribution,  while  the  other  terms  will  have  a  distribution  of  the 

product  of  two  Gaussian  variables.  Again  applying  the  central  limit 

theorem  on  the  sum,  G(f  )  can  be  approximated  as  a  Gaussian  variable 

m 

with  mean 


“  f  •?  <w-f»> 
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2N  Tm 


-  N  WT  m  (1  -  f  /W) 
o  m 


and  variance 
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2  2  4 

2N  T  m  ' 


-  2N  2T3W  (1  +®“]2  [1  -  f  /V]  (40) 

o  /  n 

If  the  terras  for  b(t)  are  taken  into  account,  both  the  mean  and 

I 

variance  of  F  (f  )  will  be  twice  what  was  calculated  to  Equations  39 

VI 

and  40.  Therefore,  the  distribution  of  F  (f  )  is  Gaussian  with  mean 


Uf  -  2M  VT  a  [1  -  f^/W) 
» 


and  variance 


°2f  ’  4n02t3w  11  *  ?“j2  (1  *  *  /w) 


Since  the  variance  in  Equation  32  must  now  be  normalised  as 


in  the  case  of  modulation,  a  similar  procedure  is  used  to  obtain 
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a  value  for  Nq  .  Referring  to  Equation  20,  the  following  equation  is 
obtained, 


4N  Vv  (1  -  f/W)  (1  +  m2/2)2  -  (1  -  f/W) 

o 


3  2  2 

4T  W(1  +  m  /2) 


Therefore,  the  distribution  of  a  sample  of  a  frequency  f 

0 

will  be  Gaussian  with  mean 


mVwf  11  -  tJV) 

(1  +  m2/2) 


and  variance 


(1  *  f/W) 


CHAPTER  IV 


EXPERIMENTAL  ANALYSIS 

The  Information  given  by  Equations  44  and  45  is  sufficient 
for  finding  the  probability  of  detection  of  the  modulation  as  a 
function  of  the  modulation  factor.  However,  to  verify  these 
equations,  they  must  be  modified  so  that  the  experimental  results 
can  be  compared  with  the  theoretical  lesulcs.  The  modification 
oi  these  equations  will  depend  on  what  equipment  Is  used  to  find 
the  experimental  values  and  this  cun  be  seen  best  by  explaining 
the  experimental  setup.  A  block  diagram  of  the  system  used  is 
shown  in  Figure  4.1. 

A  modulating  freque<  uy  of  50  Hz  was  used  to  modulate 
broadband  noise.  The  modulated  noise  was  then  passed  through  a  2*4 
kHz  octave  band  filter  and  amplified.  To  obtain  the  square  of  the 
envelope,  the  modulated  noise  was  rectified,  squared,  and  passed 
through  a  low  pass  filter  The  resulting  signal  was  then  converted 
into  the  frequency  doe*  In  by  use  of  a  Quan-Tech  Wave  Analyzer.  A 
one-cycle  bandwidth  filter  In  the  analyser  was  used  In  keeping  with 
the  equations  derived  in  Chapter  111.  The  analyzer  was  set  on  the 
frequency  of  interest,  50  Hi,  and  the  amplitude  of  the  output  of 
the  analyser  at  that  frequency  was  recorded  on  an  X-Y  plotter  for 
500  seconds.  The  amplitude  distribution  of  the  output  of  the 
analyzer  was  found  by  plotting  a  histogram  of  the  amplitude  at 
every  half -second  for  the  500  seconds.  The  curve  beceaM 
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sufficiently  smooth  after  using  1000  points  so  that  it  could  be 

compared  with  a  theoretical  curve. 

A  50-second  sample  of  an  X-Y  recording  is  shown  in  Figure  4.2. 

In  this  sample,  no  modulation  was  present  on  the  input  signal.  Since 

the  amplitude  of  the  output  voltage  ranged  from  i.ero  to  0.32  volts 

over  the  entire  500  seconds,  amplitude  increments  of  0.02  volts 

were  used  to  plot  the  histogram.  For  each  modulation  factor,  ten 

samples  of  data  similar  tc  that  shown  in  Figure  4.2  were  used  to 

find  the  amplitude  distribution  of  the  output  of  the  analyzer 

The  theoretical  distribution  can  be  fcund,  knowing  how  the 

spectrum  analyzer  treats  the  signal  In  this  case,  the  output  of 

the  one-cycle  filter  in  the  analyzer  is  passed  through  a  voltage 

doubler  and  a  low  pass  filter  In  effect,  this  gives  the  envelope 

of  the  signal  coming  out  of  the  filter  In  Chapter  111,  the 

output  of  the  one- cycle  filter  was  found  to  have  a  Gaussian 

distribution  with  mean  p  given  by  Equation  44  and  variance  j  1 

8  S 

given  by  Equation  45.  If  there  is  no  modulation  present,  p  is 

3 

zero  and  the  distribution  of  the  Gaussian  noise  will  be  a  Rayleigh 

distribution.  When  modulation  is  introduced,  the  mean  p  becomes 

s 

nonzero  and  the  distribution  of  the  envelope  will  be  a  function  cf 

u  Rice  has  derived  the  distribution  o 1  the  envelope  of  Gaussian 
8 

ncise  plus  a  D.C.  voltage  In  this  case,  the  D  C  voltage  is  the 
mean  of  the  Gaussian  distribution  The  distribution  of  the 
envelope  of  Gaussian  noise  with  mean  p#  and  variance  ,  according 
tc  Rice  is  the  following. 
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IS  2  E„ 

p  (E)  *•  exp  ~  - -  Iq  j  — |  [  ,  »46) 

a  2  a  i  '  7 

S  S  •  8  1 

where  I  (xj  is  the  modified  Bessei  function  of  zero  order  and  can 
o 

be  conveniently  found  in  the  Chemical  Rubber  Company  Tables  [5j. 

To  normalize  the  above  equation  with  the  distributions  found 

experimentally,  the  mean  of  the  -I  •>ve  equation  for  n8  *  0  was  set 

equal  to  the  mean  found  in  the  experimental  distribution  for  n  •  0. 

The  mean  of  Equation  46  for  ■  0  is  1.2S  o  so  that  the  value  of 

ag  can  he  calculated.  The  experimental  mean  for  the  distribution 

2 

was  found  to  be  0.118v.  so  that  o  *  0.0944  and  a  *  0.0089. 

s  s 

Knowing  this  value  and  calculating  ^  from  Equation  44  for  a  given 
value  of  m  and  normalizing  pg,  Equation  46  can  be  plotted  and  compared 
to  the  distributions  found  experimentally.  Five  different  modulation 
factors  of  ■  *  0,  0.02S,  0.05,  0.075,  and  0.1  were  used  in  the 
experiments  and  the  plots  of  the  distributions  for  each  of  these 
modulation  factors  is  shown  in  Figures  4.3  -  4.7.  The  indicated 
modulation  factors  for  the  experimental  distributions  were  found  by 
using  the  mean  of  the  distribution  to  determine  what  value  of  m 
would  give  that  mean  theoretically.  These  indicated  modulation 
factors  are  well  within  the  accuracy  of  the  equipment  used  for  the 
experiment . 

With  the  information  available  from  the  preceding 
distributions,  confidence  percentsges  can  be  calculated  for  a  given 


Figure  4.6  Analyzer  Output  Distribution  for 
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modulation  factor  The  confidence  percentage  Is  defined  as  the 
probability  of  detecting  the  modulation  for  a  given  threshold  level. 
The  threshold  level  can  be  calculated  for  a  given  probability  of 
allowing  the  signal  to  exceed  the  threshold  for  zero  modulation 
Since  the  amplitude  distribution  of  the  signal  for  zero  modulation 
was  found  to  be  a  Rayleigh  distribution,  the  thresholds  can  be 
calculated  from  the  following  equation 


where  t  is  the  threshold  and  p  is  the  probability  of  the  signal  being 
above  t  The  left-hand  side  of  Equation  47  is  mtegrable  in  closed 
form  and  becomes* 


1 


P 


«-'2/2’.2 


*  *»8  • 


Since  a ^  is  known  to  be  0  0089 ,  values  of  t  can  br  calculated  for 
given  values  of  p  and  these  values  are  shown  in  Tsble  4  . 

For  each  of  these  thresholds  s  confidence  percentsge  csn  be 
calculated  by  integrating  Equation  46  from  t  to  •  This  wii,  give 
the  probability  that  the  modulated  signal  is  above  the  threshold  t 
Since  the  expression  for  the  distribution  cannot  be  integrated 
readily,  the  Integration  was  done  numerically  by  the  use  of  Simpson's 
rule.  Figures  48  -  4.11  show  the  result  of  this  integration  as 
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TABLE  4.1 


PROBABILITY  OF  THE  SIGNAL  BEING  ABOVE  THE  THRESHOLD 


Threshold  t 


Probability  of 
Signal  Being  Above. t 


0.28 


0.01 


0.22 


0.05 


0.20 


0  10 


0.16 


0  25 


lonfidence  Percentages 


Figure  4.11  Confidence  Percentages  for  t  ■  0.16 


a  function  of  modu  c.tion  factor.  The  experimental  results  were 
found  from  the  histograms  by  counting  the  number  of  samples  aDcve 
the  threshold  and  dividing  by  1000  This  gives  an  approximate  value 
for  the  area  of  the  curve  above  the  threshold.  Prom  these  plots,  the 
percentage  of  the  time  the  signal  is  above  a  given  threshold  can 
be  found  as  a  function  of  the  modulation  factor. 

It  must  be  remembered  that  the  experimental  results  in  this 
study  are  obtained  for  noise  of  a  particular  bandwidth  of  2  kHz 
It  can  be  seen  from  Equation  44  that  the  mean  value  of  the  spike  at 
the  modulating  frequency  is  proportional  to  the  square  not  of  the 
noise  bandwidth.  Suppose  the  noise  was  of  bandwidth  20  kHz. 

This  would  increase  the  seer,  value  of  the  spike  by  a  factor  cf 
3  lb  and  would  make  the  ccnfldence  percentages  even  higher  sin:* 
the  output  of  the  analyzer  would  be  above  the  threshold  more  of  the 
time  Likewise,  the  confidence  percentage  for  a  given  modulation 
factor  would  be  lesa  if  the  bandwidth  of  the  notes  Is  laaa  than  the 
2  kHz  used  in  the  experiment.  Tha  contidenoe  percantages  for 
different  bendvidths  can  be  found  by  rsplotting  Equation  4b  with  the 
new  value  for  the  mean  and  by  finding  the  area  under  tha  curve 


aoove  the  threehold  value 


CHAPTER  V 


SUMMARY 


Results  and  Concluslom 


An  important  result  of  this  study  Is  found  in  the  curves  of 

Chapter  IV  showing  the  percentage  of  time  the  frequency  spectrum 

of  the  envelope  of  modulated  Gaussian  noise  is  above  a  certain 

threshold  at  the  modulating  frequency.  It  can  be  seen  from  these 

curves  that  10  percent  modulation  will  almost  always  be  above  a 

threshold  for  which  a  signal  with  no  modulation  would  only  be  above 

one  percent  of  the  time.  This  means  that  the  spike  at  the  modulating 

frequenc>  In  the  frequency  spectrum  will  be  easily  detectable 

anove  tht  adjacent  frequencies  for  any  noise  of  bandwidth  2  kHz 

that  is  modulated  at  least  10  percent  by  a  sine-wave.  Of  course, 

lowering  the  threshold  will  increase  the  chance  of  detecting  the 

spike  but  that  also  will  Increase  the  false  alarm  rate.  That  is, 

there  will  be  a  greacer  chance  of  the  signal  being  above  the 

threshold  when  there  is  no  modulation  present. 

A  more  general  result  that  has  been  found  is  that  the 

distribution  of  the  frequency  spectrum  for  the  envelope  of  modulated 

Gaussian  noise  is  a  Gaussian  distribution  with  mean  p  given  by 

8 

2 

Kcuaticn  44  and  a  given  by  Equation  45.  In  this  study,  the 

8 

spectrum  analyzer  used  gives  the  envelope  of  this  frequency 
spectrum  but  this  is  not  always  the  case.  For  example,  several 
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spectrum  analyzers  use  a  square  law  detector  at  the  output  to 
give  the  energy  spectrum.  This  would  give  an  output  entirely 
different  trcm  the  one  that  was  obtained  in  this  study.  In  any 
:a se..  knowledge  of  hew  the  spectrum  analyzer  treats  the  signal  Is 
essential  m  determining  the  amplitude  distribution  of  the  output 
of  the  analyzer  and  this  distribution  must  be  found  so  that  the 
theoretical  and  experimental  results  can  be  compared. 

5.2  Areas  for  Further  Study 

One  area  for  further  study  would  be  to  find  the  probability 
of  detection  of  modulation  by  using  the  amplitude  distribution 
of  the  modulated  Gaussian  noise  rather  than  the  spectral 
distribution  as  used  in  this  study.  The  results  could  then  be 
compared  with  the  probabilities  of  detection  found  in  this  study. 

Another  area  for  study  would  be  to  consider  several  other 
types  of  modulation  beside  jine  wave  modulation.  Amplitude 
modulation  of  a  noise  carrier  from  a  radiated  noise  source  may 
not  aiways  be  sine  wave  modulated.  The  modulation  might  be  closer 
to  square  wave  modulation  or  even  triangular  wave  modulation.  These 
types  of  modulations  should  be  studied  and  compared  with  the  results 
of  this  study  which  Involves  only  sine  wave  modulation. 
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